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Abstract 



Q^' We combine, within the framework of the Standard Model, the results of 

o,' Higgs search experiments with the information coming from accurate the- 

O , oretical calculation and precision measurements to provide a probability 

density function for the Higgs mass, from which all numbers of interest can 

be derived. The expected value is 170 GeV, with an expectation uncertainty, 
S^ . quantified by the standard deviation of the distribution, of about 80 GeV. 

H \ The median of the distribution is 150 GeV, while 75 % of the probability 

is concentrated in the region Mh < 200 GeV. The 95 % probability upper 

limit comes out to be around 300 GeV. 



1 Introduction 

Presently, one of the main interests in High Energy physics is the search for 
evidence of the Higgs boson and the determination of its mass. Although all 
direct searches have been unsuccessful till now, the self consistency of the 
Standard Model (SM) in the electroweak sector [|^] makes physicists highly 
confident about the hypothesis that the Higgs boson exists, and, most likely, 
it has effective properties close to those expected from the minimal Standard 
Model. 

In this paper we use the accurate theoretical predictions for the effective 
mixing parameter, sin"^ O^Fr = s^rr Q, and the W boson mass, Mw ^, 
together with available experimental information, including also the results 
of direct search experiments carried out at LEP, to infer the value of Higgs 
mass, Mh- Clearly, the unavoidable status of uncertainty on the value 
of each of the experimental parameters, as well as on the accuracy of the 
calculations, allows only a probabilistic inference to be made. As a result 
we provide a probability density function (p.d.f.) for the mass of the Higgs 
boson 

firuH I "data","SM" ) = f{mH \ dir. & ind.) 

conditioned by the experimental data from direct searches [dir.) and pre- 
cision measurements {ind.) under the assumption of validity of the SM. 
Prom this function we make a set of probabilistic statements about Mj/, 
and summarize the result in terms of convenient and conventional numbers 
(expected value, standard deviation, mode, median, etc.). 

The paper is structured in the following way: in the next section we 
recall the theoretical formulae used in the analysis; section ^ is devoted to 
a detailed description of the inferential method used; all the input quanti- 
ties which enter the analysis are introduced and commented on section ^. 
Section |5| deals with the determination of Mw and Mh using only indi- 
rect information. Section |6| presents the main result of the paper, namely 
f{mH I dir. & ind.). Finally we draw some conclusions. 

2 SM formulae relating the Higgs mass to the ex- 
perimental observables 

The most convenient way to approach the problem is to make use of the 
simple parameterization proposed in Ref.[Q], in which the relations among 
the observables mostly sensible to the Higgs boson mass are summarized in 



two formulae: 






«2 


= i^lff)o + ClAi + C2A2 - C3A3 + C4^4 , 


(1) 


Mw = 


= M^ - diAi - d^AJ - d2A2 + ^3^3 - diAi . 


(2) 



{slff)o, M^, Ci and di are theoretical quantities and Ai are related to experi- 
mental observables, namely Ai = ln(MH/100 GeV), A2 = [(Aa)?,/0.0280 - 1], 
A3 = [(Mi/175GeV)2-l] and A4 = [{as{Mz) /0.118 - I], where Mt is 
the top quark mass, as{Mz) is the strong coupling constant and {Aa)h is 
the five-flavor hadronic contribution to the QED vacuum polarization at 
q^ = M|. The two theoretical quantities {s^f.ff)o and M^ are the analogues 
of the experimental ones, but obtained by the theory at the reference point 
{/\a)h = 0.0280, Mt = 175 GeV, and as{Mz) = 0.118. 

Eqs. (y) and (y) are simple analytical formulae that reproduce to very 
good accuracy the results of Refs.[||, Q. In these papers the incorporation 
of the 0{g^Mf /My^) corrections into the calculation of the effective elec- 
troweak mixing angle and M\y is presented. These new contributions are 
now implemented in the fitting codes used by the LEP Electroweak Working 
Group (EWWG). Their analysis shows that the effect of these new correc- 
tions is to lower the fitted value for the Higgs mass by about 30 GeV and to 
cause a sizable decrease in the ambiguity related to the scheme dependence 
[^. Tables 1 and 2 of Ref.[Q] report the values of the various theoretical 
quantities entering in Eqs. (||) and (||) for three different renormalization 
schemes. Using these values, Eq. (|l|) approximates the detailed calculations 
of Refs.y, ^ for 75 < Mh < 350 GeV, with the other parameters in the 
ranges around the reference values 5Mt = ±5.5 GeV, 5{IS.a)h = ±0.00065 
and 5as = ±0.05, with average absolute deviations of ~ 4 x 10~^ and maxi- 
mum absolute deviations of (1.1 — 1.3) x 10~^ depending on the scheme. Sim- 
ilarly Eq. (H) shows average absolute deviations of approximately 0.2 MeV 
and maximum absolute deviations of (0.8 — 0.9) MeV. Outside the above 
range, the deviations increase reaching (2.6 — 2.8) x 10~^ and (3.1 — 3.3) MeV 
at niH = 600 GeV. It is clear that augmenting Eqs. (jl]) and (^) with higher 
powers of the Ai quantities one can reach a better agreement in a wider 
range of Higgs mass values. However, once "a posteriori" it is verified that 
indeed /(m-jy) is concentrated in the region of validity of Eqs. (ffl) and (0) 
then the use of more complicated parameterization is not really necessary. 

Eqs. (|l|) and (0) can be used for two different purposes. The first is to 
solve them with respect to Ai and M\y, getting simultaneously Mh and 
M\Y. This value of M\y can be compared to the experimental measure- 
ments in order to check the consistency of the theory. Once this check is 



performed, both Eqs. (|T]) and (Q) can be used to infer Mh-, although the 
two determinations are not independent, due to common terms in the two 
formulae. The combination of these two results, taking into account the cor- 
relations, provides a joint distribution that, further constraint by the direct 
search results, will give us the final f{mH \ "data" , "SM" ) . 

3 Analysis method 

Probabilistic approach to infer Mh and Mw 

The quantities entering Eqs. (|l|) and (^) (hereafter called "input quantities") 
are not known exactly, and this makes the result uncertain too. It is natural 
to handle this uncertainty by probability. The numerical value of the input 
quantities, which here will be generically indicated]^ by Xi, is interpreted as 
an uncertain number (also called "random variable"). This means that each 
of them could assume an infinite number of possibilities, each characterized 
by a number /(xj), such that f{xi)dxi gives the (infinitesimal) probability 
that the "true value" Xi is in the interval dxi around Xi. In this framework 
the extraction of Mh from Eq. (|l|) gives a solution that depends on the 
uncertain values Xi, 

MH = MH{Xi,X2,...,Xn), (3) 

and therefore we need to evaluate the p.d.f. of a function of random vari- 
ables. The most general way to describe the uncertainty about the value 
of quantities Xi is given by the joint distribution /(xi, X2, ■ ■ ■ , x„). Then, a 
straightforward application of the probability calculus leads toH: 



fimn) = / f{xi,X2,...,Xn) ■ S{mH - Mh{xi,X2, . . . ,Xn)) dxidx2 ■ ■ ■ dxn , 

(4) 
where the integral is extended over the hypervolume in which Xi are defined. 
The l.h.s. of Eq. (|) actually stands for /(m^ | ind., "Eq. 1^"). Eq. (|) has 
a simple intuitive interpretationg: the (infinitesimal) probability element 
f{mH) dniH depends on "how many" elements dxidx2 ■ • • dxn contribute to 
it, each weighted by the p.d.f. calculated in the point {xi, X2, . . . , x„}. 

^Notice that, following the practice of probability theory, we indicate with capital 
letters the name of the variable and with small letters the values they may assume. 

^An alternative way of interpreting (m is to think to a Monte Carlo simulation, where 
all the input quantities enter with their distributions, and correlations are properly taken 
into account. The histogram of Mh calculated from (H) will "tend" to fijnH) for a large 
number of generated events. 



The solution of Eq. (^ is very complicated, however we can perform a 
series of approximations and use of central limit theorem to get the final 
p.d.f. without actually making explicit use of Eq. (^ and without reducing 
the accuracy of the inference. We would like to list the steps needed to 
determine fimu \ ind., "Eq. ([l|)"). 

• First, with a great degree of approximation, the quantities entering 
Eq. (|l|) are independent, or at least this condition is satisfied for the 
quantities from which the uncertainty on Mh mostly depends. Ac- 
tually, the theoretical parameters entering Eqs. (Q) and (Q) contain 
the same information evaluated in different renormalization schemes, 
and, therefore, they could all be correlated. A more careful procedure 
for handling their uncertainty could be considered. This issue will be 
discussed at the end of this paragraph, and the numerical outcomes of 
the two methods used will be compared when discussing the results. 

• Second, we make use of the central limit theorem, which makes the 
probability distribution of a linear combination of random quantities 
under well known conditions Gaussian. The importance of this theo- 
rem is that we only have to make sure that the terms dominating the 
overall uncertainty are practically Gaussian. As far as the other terms 
are concerned, the exact form of the individual distributions doesn't 
even matter, since only expected value and variance are relevant. 

• The consequences of the central limit theorem can be extended to the 
variables which do not enter linearly, if their dependence can be lin- 
earized with a reasonable degree of approximation in a range of several 
standard deviations around their expected value. This amounts to re- 
quiring these variables to have a sufficiently small variation coefficient 
(the "relative uncertainty" of the physicists' lexicon). 

• Applying this analysis to our case, we see that the solution of Eq. (|l|) 
in terms of tjih is strongly not linear. Therefore Ai is the natural 
quantity with which to express the result at an intermediate stage, 
being 

Ai-AA(E[^i],a(Ai)), (5) 

where the last notation is a shorthand for a normal distribution of 
expected value E[^i] and standard deviation a{Ai), calculated as 

E[^i] = A^{E[Xi],E[X2],...,E[Xn]) (6) 



with the derivates evaluated at the expected values. 

• Finally, the exact form of f{m}{) can be obtained from /(ai), making 
use of standard probability calculus, e.g. using Eq. (Q). 

A similar strategy can be used to get the parameters of the Gaussian which 
describes the knowledge of My/. In this case the linearization hypothesis 
is already reasonable for My/ itself and the resulting f{mw) is therefore 
normal with a high degree of approximation. 

The procedure outlined above does not take into account possible cor- 
relations among the theoretical parameters of Eqs. (^ and (pf). Estimating 
the correlation coefficients from the sample provided by Tables 1 and 2 of 
Ref.[Q] would be a rough and complicated procedure. In fact the covariance 
matrix can only take into account linear correlations, whether, in general, 
these effects could be more subtle. A more elegant and general way to handle 
this information is, then, to consider different inferences, each conditioned 
by a given set of parameters, labelled by Ri. This can be applied at any 
stage of the analysis, although it is in practice more convenient to apply it 
at the level of the inference on Ai. For each renormalization scheme Ri we 
have then: 

Ai\R^^M{nAi\Rila{Ai\Ri)). (8) 

The p.d.f. of j4i, "integrated" over the possible schemes, is then 

f{ai) = Y.f{ai\R,)-f{Ri), (9) 

i 

where f{Ri) is the probability assigned to each scheme. The calculation of 
expectation value and variance is straightforward. When there is indifference 
with respect to the renormalization schemes (i.e. f{Ri) = 1/3 Vz) we get 

E[^i] = ^Y.^[Ai\R,] (10) 

i 

a\A,) = i^a2(Ai|i?,) + i^E2[Ai|i?,]-E2[^i] (11) 

i i 

= i^a2(^i|i?,)+4, (12) 

i 

where ge indicates the standard deviation calculated from the dispersion of 
the expected values. The p.d.f. (0) is in general not Gaussian, since it comes 
from a linear combination of p.d.f. 's, and not from a linear combination of 
variables (i.e. the central limit theorem does not apply). Nevertheless, in our 
case the Gaussian approximation will be valid, as will be discussed below. 



Double inference on Mh and combination of the results 

The method described in the previous section is apphed to each of the 
Eqs. (^ and (^), obtaining two inferences on ^i, the first (indicated by 
j4|) depending on the effective electroweak mixing parameter and the sec- 
ond {A^) on the W mass. The second equation leads to two solutions and 
the largest value has been considered, because of the agreement with the ^f 
and also because the smaller solution leads to a mass well below the range 
firmly excluded by present observations. 

The two uncertain values A\ and A^ are not independent, due to the 
fact that some of the input quantities appear in both relations. This means 
that we have to consider the joint distribution f{a\,ai). Because each 
variable is individually Gaussian, the joint distribution is described by a 
two-dimensional normal, with a correlation coefficient p{A\,A^) calculated 
from the covariance between A\ and A^: 



<r(Af).<r(Ar 



Again using linearization around the expected values, one finds easily that 
the covariance is given by 






dXi dXi 

Eflf-W))-flt^-(^.)). (13) 



^ \dXi ' '7 \dXi 
where this last formulation is very convenient for practical purposes, as we 



will see below. Eq.(13) does not take into account the correlations among the 
various theoretical coefficients. However, numerically they are completely 
negligible with respect to the ones due to A2, A^, and ^4 and therefore the 
use of Eq. ( [T3D is well justified. 

The presence of the correlation term prevents the two results from being 
combined with the usual formula of the average weighted with the inverse 
of the variance. There are several possibilities for taking correlation into 
account, either working directly with p.d.f.'s, or assuming that the final 
result is also normally distributed and evaluating the two parameters of the 
distribution. Obviously the conclusions will not depend on the procedure 
if the normality assumption is correct, as it is in this case. The way that 
seems to us the most intuitive relies on the fact that p{Ai,Ai) is positive, 
as we will see, and thus the correlation between the two results is equivalent 



to that introduced by an uncertainty on a common offset (see, e.g., |^). 
Therefore the variances of Af and Af may be considered as being formed of 
two parts: one of these parts, indicated by o"^, is common to both variances; 
while the other is individual. The common part is given by covariance, i.e. 
(7^ = Cov(Al,Af). The individual contribution to each variance is then 
evaluated subtracting o"^. This procedure allows the expected value E[^i] 
to be evaluated as the average of E[Af ] and E[^]"], weighted with the inverse 
of the individual variance. The variance a'^{Ai) will be, finally, the sum of 
the "variance of the weighted average", plus the common variance.^ We 
have, then: 
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Finally, the Gaussian result on A\ is transformed into the p.d.f. of Mu 
using probability calculus. For example, one can make use of Eq. (Q), 
and the result is straightforward, remembering that bijxiu — 100 expai) = 
roll (5(ai— ln(m,j:f/100)). Then it is possible to evaluate expected value, stan- 
dard deviation, mode (Mh) and median (Mfj) of Mh (see, e.g., Q for the 
properties of the so called lognormal distribution). The results (expressed 
in GeV) are: 



f{mH 1 ind.) = 


1 


1 


(ln(m^/100)-E[^i])2" 
2a2(Ai) 


(14) 


V2 7ro-(Ai) rriH 


E[Mh] = 


= 100 exp 


E[A,] + -a\A,) 


(15) 


a{MH) = 


1 
= 100 (exp[2E[Ai] + 2fT2(Ai)] - exp [2E[Ai] + ^^(^i)]) "(16) 


Mh = 


= 100 exp 


E[Ai] - a\Ai)] (17) 


M^o = 


= 100 exp 


E[^i]] 






(18) 



Notice that the value of these position and dispersion parameters of /(ttt.//) 
is, in general, not simply the back transformation of those of /(ai). 

^An alternative way, which still avoids working with p.d.f. 's, is described in M. The 
two procedures yield identical results. 



Including the constraint from direct search 

The knowledge about the value of Mh is modified further by the non- 
observation of the Higgs boson up to the highest LEP energies. To un- 
derstand how f{mH I ind.) changes when it is further conditioned by the 
negative direct experimental result, let us consider a search for Higgs pro- 
duction in association with a particle of negligible width in an ideal situation 
( "infinite" luminosity, perfect efficiency, no background) whose outcome was 
no candidate. Consequently all mass values below a sharp kinematical limit 
Mk are excluded. This implies that: a) fimn) must vanish below Mk 
(otherwise one would have observed the particle); b) above Mk the relative 
probabilities cannot change, because there is no sensitivity in this region, 
and then the experimental results cannot give information over there. For 
example, if Mk is 90 GeV, then /(200GeV)//(100GeV) must remain con- 
stant before and after the new piece of information is included. In this ideal 
case we have then 

r mn < Mk 

f{mH I dir. & ind.) = I f{n,H\ind) ^ > ^ (19) 

[ /^7 f(mH\ind.)dmH " - ^ ' 

where the integral at denominator is just a normalization coefficient. 

More formally, this result can be obtained making explicit use of the 
Bayes' theorem. Applied to our problem, the theorem can be expressed as 
follows (apart from a normalization constant): 

firriH I dir. & ind.) oc f{dir. \ mu) ■ f{mH \ ind.) , (20) 

where f{dir \ mn) is the so called likelihood, which has the role of updating 
the p.d.f. once the new piece of information is included in the inference. In 
the idealized example we are considering now, f{dir \ niH) can be expressed 
in terms of the probability of observing zero candidates in an experiment 
sensitive up to a Mk mass for a given value mn, or 

f{dir. I mu) = /("zero cand." | mu) = < ^ m^ > M^ ^'^^^ 

In fact, we would expect an "infinite" number of events if Mh were below 
the kinematical limit. Therefore the probability of observing nothing should 
be zero. Instead, for Mh above Mk, the condition of vanishing production 
cross section and no background can only yield no candidates. 

In real life situations the transition between values which are impossible 
to those which are possible is not so sharp. Because of physical reasons 



(such as threshold effects and background) and experimental reasons (such 
as luminosity (£) and efficiency (e)) we cannot be really sure about excluding 
values just below Mk, nevertheless very small values of the mass are ruled 
out. In the case of Higgs production at LEP the dominant mode is the 
Bjorken process e^e~ -^ H + Z° . This reaction does not have a sharp 
kinematical limit at ^/s — Mz (minus a negligible kinetic energy), due to 
the large total width of the Z°. The effective kinematical limit {Mk^..) 
depends, then, on the available integrated luminosity and could reach up 
to the order of ~ ^/s — Mz + 0(10 GeV) for very high luminosity. This is 
clear from figure |l|a, where the cross section e"'"e~ -^ H + Z°, with the Z 
decaying in all possible channels, is plotted as a function of the Higgs boson 
mass for 172, 183 and 189 GeV cm. energy, with the vertical lines showing 
^/s — Mz for the three cases. 



In this non ideal situation we expect the step function of Eq. (21) to be 
replaced by a smooth curve which goes to zero for low masses and to 1 for 
Mh — > Mx^ff, where the experimental sensitivity is lost. In this case dir. in 
Eq. (|20| ) stands, in principle, for all possible experimental observables and 
the function f{dir. \ mn) should be provided by the experiments. However 
the LEP collaboration results on the Higgs mass searches are usually pre- 
sented in terms of confidence level (C.L.) [|lO|. As discussed in section |6| this 
quantity does not have a simple connection to f{dir. \ mn)- Given this sit- 
uation we decided to model the likelihood in a way which seems compatible 
with the physics case. 

First, all possible experimental observables can, in practice, be replaced 
by suitable combinations which depend on the Higgs mass. The simplest 
of these possible "summaries" of the data is the number of observed can- 
didate events, which we will indicate by Uots- The number of candidate 
events expected to be observed, on the other hand, is given by the sum of 
the Higgs events, indicated by Xs^ttih), and the expected number of back- 
ground events, A^, assumed here to be well known (see e.g. Ref. P] for the 
natural extension when A^ is uncertain too). The mass dependence of the 
former is due to the mass dependence of cross section, branching ratio (b.r.) 
and efficiency, and so it depends on the decay channel investigated. For 
simplicity, we discuss the case of a likelihood obtained considering the total 
number of observed candidate events in a single channel. This is given by 

„. , ^ , e-(^s(m^)+AB).(A5(,nH) + AB)"°*^ ,^^, 

f[nobs\mH,XB) = j , (22) 

since Uobs is expected to be described by a Poisson distribution with pa- 
rameter A = Xsi''TT-H) + A^. In order to compare and combine the updating 
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Figure 1: a) cross-section e+e — > H + Z° as a function of ttih for 
^/s = 172, 183, 189 GeV. The vertical lines are at mn = ^/s - Mz- b) 7^ 
vs. ruH for Uots = As = for ^s = 172, 183, 189 GeV, with e = 30% and 
£ = 10, 55, 180 pb^^ per experiment, respectively. Odd lines are for a single 
experiment, even ones represent the combination of four. 
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power provided by each piece of information easily, it is convenient to re- 
place the likelihood by a function, TZ, that goes to 1 where the experimental 
sensitivity is lost [^I|. This function can be seen as the counterpart, in the 
case of a real experiment, of the step function of Eq. (^T]). Because con- 
stant factors do not play any role in the Bayes' theorem we can divide the 
likelihood by its value calculated for very large Higgs mass values, where no 
signal is expected^, i.e. Mh — > "oo", or A^ — > 0. Clearly, this operation 
makes sense only if the likelihood is different from for Mh -^ "oo" . This 
condition is satisfied for any A^ in case Uobs / 0, but when Uobs = only for 
Xb = 0. The case of Uobs 7^ with A^ = leads to a clear discovery, i.e. the 
likelihood will assume positive values only below Mk^,, and there is no need 
anymore to build the TZ function with the desired asymptotic properties. 

IZ, as a mathematical function of niH-, with Uobs and Xb acting as param- 
eters, can be seen as a kind of shape distortion function of the p.d.f. intro- 
duced by the new data. As long as TZ{mH) is 1, the shape (and therefore the 
relative probabilities in that region) remains unchanged, while in the limit 
TZ{mH) -^ the p.d.f. vanishes. One should notice that TZ{mH) can also 
assume values larger than 1 in the region of sensitivity, corresponding to a 
number of observed candidate events larger than the expected background. 
In this case the p.d.f. will be stretched below the effective kinematical limit 
and this might even prompt a claim for a discovery if 7^ becomes sufficiently 
large for the probability of Mb in that region to get very close to 1. 

Applying this formalism to our case and in the realistic situation of non 
vanishing expected background we get 



-(As("xh)) . (X^mH) + XbT°'" 



n{mH]nobs, Xb) = ni^j^ , (As / 0, if nobs + 0). 

\^ 

(23) 

Instead, when As = and n^bs = one can take the limit of the above 
formula obtaining TZimn) = e"^''™'^^ Examples of this function are 
shown in figure |l]b) in case of no events candidates and zero background 
iP'obs = As = 0), for y^ = 172, 183 and 189 GeV, considering a single LEP 
experiment (odd numbers) and the combination of all experiments (even 
numbers). The calculations have been done assuming a nominal integrated 
luminosity per experiment of 10, 55 and 180 pb~^ for the three energies, 
and an average and constant detection efficiency of 30%. Fig.g, instead, 
illustrates six different possible scenarios. In the figure we consider a search 
at i/i = 183 GeV with £ = 55pb~^ by a single experiment that looks for 

*In the statistics lexicon this function is the Bayes factor between the generic mass 

niH and A^h = "oo" . 
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Figure 2: TZ vs. mn for two search channels at ^/s = 183 GeV with C = 
55pb~^. The lines a-c represent a channel with b.r. x e = 1.7% while d-f 
correspond to a search with b.r. x e = 32%. The cases considered are: Uobs = 
2, Xb = 0.2 (a, d); Uobs = 3, Ab = 3 (b, e); Uobs = 0, \b = 1-5 (c, f). 



two Higgs decay channels with different branching ratio. For each channel 
we plot 3 different situations of expected background and observed number 
of events. 

As already pointed out, f{dir. \ ttih)-, and hence Tl{mH] dir.), might be 
more complicated than the simplified likelihood used here, and it can only 
be provided by experiments. Either of the above functions would be the 
most unbiased way of reporting the experimental result and it would allow 
several pieces of experimental information to be easily combined. In fact, 
when individual experiments or decay channels are independent the overall 
likelihood is simply the product of the individual likelihoods and therefore 



nimn; "all data") = UiUiimH; "data i") , 



(24) 



and this can be used in Eq. ( |20| ) to get the distribution of Mh which takes 
into account all available data. 
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4 Input quantities entering the indirect determi- 
nation 

In this section we discuss the experimental and theoretical inputs used in 
our analysis. 

Hadronic contribution to QED vacuum polarization 

The QED coupling at the Z boson mass scale plays an important role in 
the prediction of mu- This fact has always stimulated a lot of activity on 
the exact determination of (Aa)/^. The most phenomenological analyses 
of it rely on the use of all the available experimental data on the hadron 
production in e"*" e~ annihilation and on perturbative QCD (pQCD) for the 
high energy tail [E > 40 GeV) of the dispersion integral. The reference 
value in this approach is ||l2|] 



(AQ)f -^ = 0.02804 ± 0.00065 . (25) 

In the recent past the hadronic contribution to the vacuum polarization 
has been the subject of several new investigations that advocate the use of 
pQCD down to energy scale of the order of 1 GeV [^, 14 1. In this more 
theory driven path, the various analyses differ on the energy value at which 
to start applying pQCD and on the amount of theoretical inputs used to 
evaluate the experimental data in the regions, like, for example, the thresh- 
old for the charmed mesons, where pQCD is not applicable. The common 
characteristic of these works with respect to the most phenomenological ones 
is to obtain a smaller central value for (Aa)h with a reduced uncertainty. 



The most stringent result of these theory oriented analyses is [14| 



(Aa)^^ = 0.02763 ± 0.00016 (26) 

that we will use in the sequel as reference value for this kind of approach. 

At the moment there is no definite argument for choosing one or other 
of the two approaches. The results are absolutely compatible to each other. 
However, the numerical difference between central values and uncertainties 
is such that it prevents an easy estimation of the effect of choosing one value 
instead of the other. For these reasons we decided to present our results for 
the values of (Aa)/j given by (Aa)^-^ and (Aa)^^ separately. 
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Top quark mass 

The value of the top quark used in our analysis is the combination of the 
experimental direct measurements reported by CDF [15| and DO |16]: 



Mt = 174.2 ± 4.8 GeV. (27) 

The Mt value obtained in the global fit of the EWWG ||5|, that uses as 
experimental value input Mj = 173.8 it 5.0 GeV, is actually a little bit 
smaller, i.e. Mt = 171.1 ± 4.9 GeV. The principal cause for this smaller 
value is connected with the remnant of the famous Rb "anomaly" . In our 
analysis we assume the validity of the SM and therefore we prefer to use the 



experimental result of Eq. (27) 



QCD strong coupling constant 

Among the various input quantities, the strong coupling constant at the Mz 
scale is the least important. In fact, QCD effects appear in the theoretical 



calculations of s^ff and Mw only at the two loop level. We use the world 



average [17| 



a^(Mz) = 0.119 ±0.002 (28) 



Effective mixing parameter sin^ T, 



t 

eff 



The effective Weinberg angle is the quantity that has the greatest sensi- 
tivity to the Higgs. Therefore its precise value is very important in deter- 
mining f{mH)- There is overall good agreement among all the measure- 
s' - 



ments although the two most precise ones, i.e. ^lr from SLD (s^ 



0.23109 ± 0.00029) and ^p'^ from LEP (s^^^ = 0.23225 ± 0.00038), are still 
about two and a half standard deviations apart. However, the continual 
raising of the SLD value during the recent years together with some reduc- 
tion of the ^p'g result has significantly improved the agreement between the 
average LEP and SLD determinations. In this situation we do not see any 
particular reason either for excluding the SLD values or for attributing to 



it a smaller weight |18]. Therefore we employ in our analysis the combined 
LEP+SLD average|] 

sin^ 6l^Jjf = 0.23157 ± 0.00018. (29) 
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W boson mass 

The present experimental information on Mw comes from the invariant mass 
of its decay products (LEP and Tevatron), from the threshold behaviour 
of the production cross section (LEP) and from the electroweak coupling 
constant in neutrino scattering (NuTeV,CCFR). The first two measurements 
can be considered a kind of "direct" determinations of the mass, in the 
sense that they are sensitive directly to it and not to a combination of other 
parameters of the SM. The combination of CDF fl^, DO pO] and LEP values 
[^ (including also the old UA2 measurement ||2l|]), gives H]: 

M^ = 80.39 ± 0.06 GeV. (30) 

The result of the deep inelastic scattering experiments can be reported using 
the quantity sin^ 0^„ = 1 — My^^/M^- In this case M\y can be extracted in 
terms of the very precise Mz value plus top quark and Higgs boson mass 



corrections [22|: 



where M^r indicates the result at the reference values {Mt = 175 GeV and 
Mh = 150 Ge\): 

M^ = 80.25 ±0.11 GeV. 

In order to make use of all available information we proceed in the following 
way. We evaluate A^ from M^ and M^^ separately using Eq. (|2|) (in the 
M^'^ case the Higgs and top dependence can be accounted for by redefining 
the theoretical coefficients di, d^, and M^). Once the compatibility of the 
two results has been established, we are allowed to combine directly the Mw 
values weighting them with the inverse of the variance. We obtain 

— ^ = 80.36 ± 0.05 + 0.0023 - 0.0057^i + O.O5IA3 (32) 

GeV 

which is the value employed in the analysis. Again, the Higgs and top de- 
pendence is taken into account by slightly modifying the relevant coefficients 
in Eq. (1). 

Theoretical coefRcients 

The various coefficients entering Eqs. (||) and (g) are not known exactly due 
to truncation of the perturbative series. This uncertainty is usually esti- 
mated comparing the results of different schemes of calculation that contain 
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all the available theoretical information at a given order of accuracy. Then 
the simplest procedure is to evaluate the best value and standard deviation 
associated to the uncertainty of each of the coefficients from the average and 
standard deviation of the values given in Ref. [Q] (when all renormalization 
schemes yield the same numerical results the standard deviation is that due 
to the rounding, i.e. unit of the least significant digit divided by vl2)- They 
are indicated in tables [l|, y and y and considered independent in the uncer- 
tainty propagation. Let us comment on the meaning and the use of averages 
and standard deviations for the coefficients. Taking as an example ci, we 
get (in units 10~^) E[ci] = 5.23 and cr{ci) = 0.04, obtained by the following 
numbersQ: 5.23, 5.19 and 5.26. This does not imply that necessarily one 
has to believe that 5.23 is really more preferred than the others, as a Gaus- 
sian distribution centered in 5.23 with standard deviation 0.04 would imply. 
One could imagine a uniform distribution ranging between 5.16 and 5.30; or 
a triangular distribution centered in 5.23 and going to zero at 5.13 and 5.23; 
or any other distribution having mean 5.23 and sigma 0.04. The final re- 
sult, relying on the central limit theorem, which, for the relative sizes of the 
standard deviations of interests ensures a fast convergence, will not depend 
on the shape of the particular distribution (they could also be different for 
different coefficients). 

It should be noticed that the values presented in Ref.|Q] do not cover 
uncertainties associated to QCD contribution in electroweak corrections. 
The dominant part of it is included in 5qcd, the relevant correction in the 
electroweak parameter Ap. The uncertainty in 6qcd will reflect itself in a 
correlated way in the various theoretical coefficients. 

To judge the effect of possible correlations in the values given in Ref.0] 
we use the method outlined at the end of Sect. 2. The more rigorous results 
derived with this procedure are practically identical to those obtained using 
average values and standard deviations of each coefficient. This is shown 
in table pi and El where the comparison of the two methods is presented. 
In the combined final Ai result we report an additional digit to test the 
accuracy due to rounding. Also the final shape of the p.d.f. of Ai obtainable 
from Eq. (^ is Gaussian with a good degree of approximation, since it is 
the average of three Gaussians (each of which is justified by the central 
limit theorem) and the closeness of their centers is much smaller than their 
widths. Given this situation, we present our result as a function of average 
coefficients and of their standard deviations, because this method shows the 
sensitivity of Ai to the various parameters in a clear way. 
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5 Results from the indirect determination 

The determinations of Mw is presented in table [l|. The two values reported 
are obtained for (Aq)^ -^ (Mw = 80.375 ± 0.027 GeV) and (Aa)^^ {Mw = 
80.366 lb 0.025 GeV). The two results are consistent and both are well in 



agreement with the experimental determinations given in (pC|) and (31). 
The uncertainty on the indirect Mw determination is still a factor ^ 2 
better than the present direct experimental result. The table also shows 
the contribution to the total uncertainty of each input quantity, with the 
sign of the derivative calculated in the reference point. This information 
allows the result to be corrected if any input quantity slightly changes in 
expected value or standard deviation. Through the entries in table || we 
can estimate the shift in the predicted central value due to unknown QCD 
effects in electroweak corrections. Indeed a variation in 6qcd introduces 
a shift on the calculated s^rr and Mw of Ss^rr Ki —1 • 10~'^ d{6QCD) M^ 
and 6Mw « 2 • 10'^ 6 {5qcd) M^ (GeV). For Mt = 175 GeV, 6{6qcd) has 
been estimated ss 5.2 • 10"'^ |2^]. This induces shifts in the values of sLr 
and AI\Y that amount to — 1.8 • 10~^ and 3.1 MeV respectively. Using table 
m one finds an additional uncertainty in the predicted Mw central value, 
6Mw ~ 1 MeV. 

The determination of Mh from the effective mixing angle and Mw sep- 
arately is presented in table § {{Aa)h = (Aa)^"^) and table ^ {{Aa)h = 
{Aa)f^^). All values are given in TeV to reduce the number of digits to 
the significant ones. The table also shows the combined determination. The 
values of ^f and Af are in agreement within uncertainty. However, the Mw 
determination is much less precise and the effect of combining it with Af 
has almost a negligible impact on the determination of the Higgs mass from 
Sg r r, also because of the correlation between them. In case of slight variation 
of the central values of the input quantities the Ai result can be corrected 
using the information provided in the tables. However the same procedure 
cannot be applied in case of changes of the standard deviations because Ai 
is obtained through a combination where the inverse of variances enters. For 
the same reason input quantities with large uncertainty are dumped in the 
combination and therefore they give a small contribution to the total Ai 
uncertainty. In the various tables, combination results are indicated by "*" 
below the relevant column. 

Among the various observables whose theoretical prediction depend upon 
Mh, given the present values of the A2-A4 quantities, sLr is by far the 
most effective in constraining Mh- Any other, like e.g. Mw or the leptonic 
width, has a very modest weight in a combined analysis. This fact justifies 
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Figure 3: Probability distribution function for the Higgs mass from precision 
measurements. The solid line is obtained employing {Aa)h = 0.02804 it 
0.00065, the dashed one with {Aa)h = 0.02763 ± 0.00016 



our choice of considering only one observable, M\y, in addition to s^pf. This 
situation will not change in the near future. In fact, a W mass as effective 
in the Mfj indirect determination as the present s^rr requires not only a 
very precise M]y result {a{Mw) < 25 MeV) but also a reduction in Mt 
uncertainty {a{Mt) < 2.5 GeV), as already pointed out in Ref.|24]. 

Figure ^ presents the p.d.f. of Mh obtained using only the indirect in- 
formation. The comparison of the two curves shows that the use of a higher 
central value for {Aa)h (i.e. (Aa)^ ) tends to concentrate more the prob- 
ability towards smaller values of Mh- This can be understood from the 
negative derivative of Mh with respect to {Aa)h, shown in tables |2| and ^. 
Indeed, in this case the median of the distribution is Mfj = 0.10 TeV while 
the analysis performed employing (Aa)^ gives for the same quantity a 
result ~ 0.3 TeV higher, which is still less than half a standard deviation 
of the distribution. It is interesting to note that although the Ai expected 
values and standard deviations in table y and ^ are different, they actually 
give very close 95% probability upper limit, Mf^. Similarly, while (Aa)^^ 
has an uncertainty that is approximately 4 times smaller than (Aq)|^ , the 
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standard deviations of the two Mh p.d.f.'s are very close too. Finally, we 
note that a variation in 5qcd of the order of magnitude discussed above 
increases Mf^ by « 10 GeV. 

6 Results including the direct search 

An extensive program to look for evidence of Higgs production in e^e~ 
collision has been pursued at LEP during the last decade. Presently results 
for Higgs searches by all four LEP collaborations are available for center of 



mass energies up to ^/s = 183 GeV |25, 26, 27, 28 1. The negative outcome 



of these searches has been reported as a combined 89.8 GeV "95 % C.L." 
lower bound p^, 50]. Unfortunately this result has no simple probabilistic 



interpretation regarding the Higgs mass [31|. The operational definition of 
the limit is expressed in terms of a test-statistic, X, based on the number of 
selected events and their distribution in a variable that discriminates signal 
from background, whose value measured in the data, Xobs, is compared 
to that obtained on the basis of a large number of "simulated gedanken 
experiments" |1C], or 

CLs+bjmH) P{Xs+b{mH) < ^obs) ,„„^ 

CLb{mH) P{Xb{mH) < Xobs) 

The C.L. for the signal + background hypothesis, CLg+b, is defined as the 
probability that the test-statistic is less than or equal to Xobs, where the p.d.f 
of Xg^b is obtained by the Monte Carlo generation of experiments in which 
a signal with mass mn is introduced in addition to the background. The 
corresponding C.L. for the signal, CLg, is then obtained by normalizing 
CLgj^b to CLb, the C.L. for the background only hypothesis, where the 
p.d.f. of Xb is obtained similarly to that of Xg+b but generating experiments 
with no signal. One can realize that the C.L. curves do not correspond 
to the likelihood of the observed data as a function of the Higgs mass (i.e. 
f{dir. I niH)). Therefore, the published combined limit cannot be translated 
into a number which could be used in an unambiguous way in our analysis. 
As outlined in section y, the information from the direct search could 
easily be incorporated into the analysis if we had the likelihood or, more 
simply, the TZ function. We note here that the test-statistic Xobs of method 



A of Ref.|lC] actually corresponds to TZ. The publication of its value as a 
function of mn would be sufficient for a complete analysis to be performed. 
However, at the moment, this information is not available. 

Although we are not in a position to make a detailed analysis, we can 
estimate the effect of the direct search results on the final p.d.f. of mn by 
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employing the simplified likelihood given in Eq. ( p2D and the public values 
of observed number of events, expected backgrounds and efficiencies. With 
respect to this a few remarks are in order: i) in the function TZ the important 
effect is always given by the data at the last energy point available. In fact, 
because the final likelihood of searches at different energies is given by the 
product of the individual likelihoods at the various y/s, in the region of Mh 
values close to Mk^,, of the highest energy point, the region which we are 
interested in, only the corresponding TZ will be relevant because those of 
lower energies are already saturated to 1. In our study we thus only use 
data at ^ = 183 GeV |^, |6|, |2^, Hj. ii) Differently from the other three 



experiments, L3 uses selection criteria that depend on mu- Correspondingly 
the number of selected events and backgrounds will be functions of mn- This 
situation is not compatible with our simplified likelihood. Therefore we do 
not consider the L3 data, in) Eq. ( [2^ ) does not contain any information 
related to the distribution of the signal and the backgrounds. At the level 
at which we can perform our analysis using only public results, shape effects 
cannot be taken into account. However, the observed number of events in 
the various channels reported by the three collaborations are always either 
zero or few units. Given this situation we expect that it is the event counting 
which gives most of the constraint. 

Figure ^ presents the TZ curves for the three experiments, ALEPH (A), 
DELPHI (D) and OPAL (O). Each curve is obtained by multiplying the 7^'s 
of the individual search channels of the experiment. The overall TZ for the 
combination of three experiments is also shown labeled as LEP3 . As already 
said, we cannot use L3 data. To try to take into account the L3 contribution 
we make the assumption that the L3 results are on average similar to those 
of the other experiments. We then roughly estimate the effect of L3, simply 
raising the combined TZ of the other three experiments to the 4/3 power. 
The corresponding curve is presented in Fig. |^ marked LEP4. We note that 
the OPAL TZ curve presents a bump which is connected to a small excess of 
the observed number of events with respect to the expected background in 
the qqH channel. This bump is not particularly significant for two reasons: 
i) it is not very pronounced and therefore when the OPAL TZ is combined 
with f{mH I ind.) the probability in the corresponding Higgs mass region is 
not particularly enhanced, ii) Our analysis is based on the event counting 
only and ignores the distribution of the signal and background. When the 
latter information is also taken into account it is not unlikely that this bump 
disappears. 

According to Eq. ( [20| ) the final p.d.f. is obtained by combining the 
p.d.f. coming from precision measurements (Fig. y) with the likelihood de- 
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Figure 4: 7?. vs. m// for the search at ^/s = 183 GeV. The curves A, D, 
correspond to the single experiment. The LEP3 line represents the combination 
of the three. The curve LEP4 is obtained by raising the LEP3 one to the 4/3 
power. 



rived from the LEP data (Fig. ^, rescaled to the function 7^, or 



f{mH I dir. & ind.) 



/q°° lZ{mH) f{inH I ind.) dm 



H 



(34) 



The result is shown in figure |5| where /{mn \ dir. &i ind.) is compared to 
/{rriH I ind.) for each assumption on (Aa)/j. To judge the sensitivity of 
f{mH I dir. & ind.) to the hkehhood we have actuahy evaluated Eq. (|34|) us- 
ing two different TVs (LEP3 and LEP4). The two resulting curves are indis- 
tinguishable as could be expected given that the two 7^'s differ practically by 
a shift of ~ 1/2 GeV. To envisage a more different case, we compare the final 
distribution presented in fig. ^ with that we obtain in our simplified analysis 
using only OPAL data. The difference in the expected value amounts to 4-5 
GeV, depending on the value of (Aa)/j chosen, while standard deviations 
are the same. The "OPAL" M^ is 3-6 GeV lower than the one reported in 
table S where our final results are presented. The closeness of the various 
quantities shows that our results do not depend critically upon the details 
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Figure 5: Probability distribution functions using only indirect information (solid 
line) and employing also the experimental results from direct searches (dashed 
one). 
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of the analysis. 

As expected, the inclusion of the direct search information in the Higgs 
mass probability analysis drifts the p.d.f. towards higher values of M//, 
changing its shape such that the probability of Mh values below 95 GeV 
drops to ~ 3%. Table ^ summarizes our analysis in terms of various con- 
venient parameters of the distribution. Expected value, standard deviation, 
mode and median are not very sensitive to the values of the hadronic con- 
tribution to the vacuum polarization. Also in both cases, 75 % of the prob- 
ability is concentrated in the region Mh < 0.20 TeV. Instead the choice of 
{Aa)h affects the tail of the distribution with (Aa)^"^ producing a longer 
one. Indeed while the use of only indirect information gives same values of 
Mjj, when the direct one is also included, the Mf^ obtained using (Aa)^"' 
is ~ 0.3 TeV higher than the corresponding number for (Aq)^ , and this 
effect is even more pronounced for Mff. 

7 Conclusions 

We have presented a method that allows the Higgs mass to be constrained 
combining together the indirect information coming from precision measure- 
ments and accurate calculations, with the results of the search experiments 
currently being carried out at LEP. The method makes use of the Bayes' 
theorem which allows the p.d.f. obtained from precision measurements to be 
augmented with the direct experimental information through the likelihoods 
of the various search channels. Such likelihoods should be provided by the 
experiments, possibly in the form of the TZ function that is very convenient 
for comparing and combining the various informations and it has also an 
intuitive interpretation because of its limit to the step function of the ideal 
case. Using the simplified form of Eq. (^) and the public results concern- 
ing observed number of events, backgrounds and efficiencies of the ALEPH, 
DELPHI, OPAL experiments for ^/s = 183 GeV, we have derived a p.d.f. for 
the Higgs boson mass that includes the direct search constraint. Several pa- 
rameters of the distribution has been reported and we have verified that the 
main conclusions do not depend to any great degree upon the detailed form 
of the likelihood. 

It should be clearly stated that our results are derived under the assump- 
tion of the validity of the SM and using as input quantities those described 
in section ^ Our analysis does not apply to different frameworks, like for 
example the Minimal Supersymmetric Standard Model except the case when 
all SUSY particles are decoupled. Concerning the input quantities, our re- 
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suits depend largely upon the experimental value for s^rr we have taken 
(Eq. (29)). It seems to us reasonable to employ the combined LEP+SLD 
value although, given the less than perfect agreement between the two most 
precise determinations, the suspicion that some measurements are affected 
by not yet understood systematic errors exists. 

The last run of the LEP machine was performed at ^/s = 189 GeV. Re- 
sults on Higgs searches at this energy are still preliminary and only DELPHI 
has presented a somewhat detailed analysis of the searches in the H + Z° 



channel |32|. In the case of negative results by the other collaborations 
similar to those reported by DELPHI, we can make a rough estimation of 
the output of the search at y^ = 189 GeV by saying that the TZ will move 
5-6 GeV towards higher values of Mh- The final /{mn \ dir. & ind.) will be 
correspondingly shifted in the same direction. 



We wish to thank F. Di Lodovico, R. Faccini and G. Ganis for useful 
communications and discussions. 
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-2.0-10-6 


-9.4-10-3 


(Aa), 


0.02804 


0.00065 


+ 13 




(Aa)fc 


0.02763 


0.00016 




+3.7 


Mi/GeV 


174.2 


4.8 


+ 13 


+12 


as{Mz) 


0.119 


0.002 


-0.60 


-0.54 


«2 


0.23157 


0.00018 


-20 


-21 


Mw/GeV 


80.375 


0.027 


- 




Mw/GeY 


80.366 


0.025 




^ 



Table 1: Mw determination. 
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Xi 


E[X,] 


a{Xi) 


1^ • ^iX^) 








"„2 » 


''Mw" 


Comb. 


('lff)o 


0.231525 


0.000015 


-0.029 




-0.024 


Cl 


(5.23 


0.04)-10-^ 


-3.1-10- 


-5 


-1.0-10-3 


C2 


(9.860 


0.003)-10-3 


-8.2-10- 


-6 


+1.3-10-5 


C3 


(2.74 


0.06)-10'3 


-1.0-10- 


-3 


-6.3 -lO-'^ 


C4 


(4.47 


0.06)-10-^ 


-9.7-10- 


-5 


-8.2-10-5 


M^/GeV 


80.3813 


0.0012 




+0.021 


+4.0-10-3 


di 


0.0578 


0.0004 




-2.0-10-3 


+5.4 -lO-'^ 


fi2 


0.5177 


0.0006 




-1.5-10-5 


-2.6-10-5 


ds 


0.540 


0.003 




-4.8- 10-* 


-3.5 -lO-'^ 


d4 


0.0850 


0.0003 




-4.5 - 10-'"^ 


-8.9-10-*^ 


ds 


0.00793 


0.00012 




-1.7-10-4 


+1.2- lO-'^ 


(Aa);, 


0.02804 


0.00065 


-0.44 


-0.21 


-0.41 


Mj/GeV 


174.2 


4.8 


+0.29 


+0.47 


+0.33 


a.(Mz) 


0.119 


0.002 


-0.014 


-0.026 


-0.017 


«2 


0.23157 


0.00018 


+0.34 




+0.29 


M^/GeV 


80.39 


0.06 




-0.82 


-0.15 


M^/GeV 


80.25 


0.11 




-0.45 

* 


-0.080 
* 


^f 


0.00 


_^ 


0.63 






AY 


+0.28 






1.07 




/ 


3=+0.34 


Ai 


0.05 -^ 




0.61 


- 


Mh/TbY 


0.13 


0.08 


{Mh = 


0.07 TeV, Mf -- 


= 0.10TeV) 



Table 2: Summary of indirect information. Ai = ln(M///100GeV), {Aa)h 
(Aq)^"'. See text for the meaning of "*" . 
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Xi 


E[X,] 


a{X,) 


14^ • ^m 








%/' 


''Mw" 


Comb. 


i4ff)o 


0.231525 


0.000015 


-0.029 




-0.026 


Cl 


(5.23 


0.04)-10-^ 


-2.1-10 


-3 


-2.3-10-3 


C2 


(9.860 


0.003)-10-3 


+7.6 • 10 


-5 


+7.0-10-5 


C3 


(2.74 


0.06)-10'3 


-1.0-10 


-3 


-8.4-10-4 


C4 


(4.47 


0.06)40-^ 


-9.7-10 


-5 


-8.9-10-5 


M^/GeV 


80.3813 


0.0012 




+0.020 


+1.8-10-3 


di 


0.0578 


0.0004 




-2.8-10-3 


+5.0-10-*^ 


^2 


0.5177 


0.0006 




+1.4-10-4 


+1.1- 10-5 


ds 


0.540 


0.003 




-4.7-10-4 


-1.6-10-4 


di 


0.0850 


0.0003 




-4.3 - 10-5 


-4.1 - 10-6 


ds 


0.00793 


0.00012 




-3.5-10-4 


+3.2-10-5 


(Aa),, 


0.02763 


0.00016 


-0.11 


-0.050 


-0.10 


Mi/GeV 


174.2 


4.8 


+0.29 


+0.46 


+0.30 


a.(Mz) 


0.119 


0.002 


-0.014 


-0.025 


-0.016 


«2 


0.23157 


0.00018 


+0.34 




+0.32 


M^/GeV 


80.39 


0.06 




-0.79 


-0.071 


M^/GeV 


80.25 


0.11 




-0.43 

* 


-0.039 

* 


^f 


0.28 


_j. 


0.46 






yir 


0.41 






1.01 




/ 


3=+0.29 


^1 


0.29i -^ 




0.455 


- 


M^/TeV 


0.15 


0.07 


{Mh = 


0.11 TeV, Mf^ -- 


= 0.13TeV) 



Table 3: Like table | but for {Aa)h 



(Aa)^. 
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Xi 


E(Xi) 


<y{x,) 


14^ • <x,) 








"MS" 


"OSI" "OSII" 


{Aa)h 


0.02763 


0.00016 


-0.10 


-0.10 -0.10 


Mt/GeV 


174.2 


4.8 


+0.31 


+0.31 +0.30 


as{Mz) 


0.119 


0.002 


-0.016 


-0.016 -0.015 


2 


0.23157 


0.00018 


+0.32 


+0.32 +0.32 


M^/GeV 


80.39 


0.06 


-0.072 


-0.070 -0.070 


M^/GeV 


80.25 


0.11 


-0.039 
* 


-0.038 -0.038 
* * 


Ai 1 Ri (MS) 


0.318 


^ 


0.456 




Ai 1 i?2 (OSI) 


0.293 


-^ 




0.460 


Ai 1 Rs (OSII) 


0.262 






0.448 


a2(Ai) 


= i^^<x2(Ai|iJ,)+4 


Ai (average) 


0.29i -^ 


0.455 



Table 4: Renormalization scheme dependence. The results shown here should 
be compared with those of table 0. 
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iAa)h = 0.02804(65) 



iAa)h = 0.02763(16) 



E[MH]/TeV 
cj(MH)/TeV 
Mh /TeY 
M^°/TcV 

P{Mh < 95 GeV) 
P{Mh < 0.11 TeV) 
P{Mh < 0.13 TeV) 
P{Mh < 0.20 TeV) 

Mfi^/TeV; P{Mh < Mf) 
Mfi^/TeV; P{Mh < Mf ) 

f Ml /TeV; P{Mh < Mi) ; 
\ Ma/TeV; P{Mh > M2) ^ 

f Ms/TeV; P(Mh < M3) 
\ M4/TeV; P{Mh > M4) 



;0.95 
;0.99 

0.16 
0.16 

;0.25 
;0.25 



{ind. 



0.13 
0.08 
0.07 
0.10 

44% 
53% 
64% 
86% 

0.28 
0.43 





3.2% 
19% 

38% 
75% 

0.33 

0.48 




0.12 
0.20 



{ind. 



0.15 
0.07 
0.11 
0.13 

23% 
33% 
48% 
81% 

0.28 
0.39 





2.3% 
16% 
34% 
76% 

0.30 
0.40 




0.12 
0.20 



Table 5: Summary of the direct plus indirect information. 
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